Weak localization of holes in high-mobility heterostructures 
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Theory of weak localization is developed for two-dimensional holes in semiconductor heterostruc- 
tures. Ballistic regime of weak localization where the backscattering occurs from few impurities 
is studied with account for anisotropic momentum scattering of holes. The transition from weak 
localization to anti-localization is demonstrated for long dephasing times. For stronger dephasing 
the conductivity correction is negative at all hole densities due to non-monotonous dependence of 
the spin relaxation time on the hole wavevector. The anomalous temperature dependent correction 
to the conductivity is calculated. We show that the temperature dependence of the conductivity is 
non-monotonous at moderate hole densities. 

PACS numbers: 73.20. Fz, 73.21.Fg, 75.70.Tj, 72.25.Rb 



Weak localization (WL) is an enhancement of 
backscattering caused by interference of waves propagat- 
ing by the same path in the opposite directions. WL of 
electrons results in corrections to conductivity anoma- 
lously dependent on temperature and magnetic held. 
Weak localization is studied in a variety of disordered 
metals and semiconductors including recently discovered 
graphene and topological insulators. In these systems, 
the carrier momentum is coupled strongly with spin or 
pseudospin which results in non-parabolic energy spec- 
trum. This leads to suppression of the backscattering 
called weak antilocalization (WAL). WAL is observed 
in graphene, for review see Ref. [TJ and in topological 
insulatorspSl the corresponding theory is developed in 
Refs. ISE 

Heterostructures with hole type of conductivity also 
belong to the class of systems with nonparabolic en- 
ergy spectrum and coupled spin and orbital degrees of 
freedom. Since the two-dimensional valence band states 
are formed as a result of competition of size quantiza- 
tion and strong spin-orbit interaction present in the bulk 
semiconductor, the spectrum nonparabolicity depends on 
the two-dimensional hole density p and the quantum well 
width a. At small density and in narrow structures the 
holes behave as ordinary electrons, while at pa 2 > 1 the 
nonparabolicity effects are important and spin is strongly 
coupled with the momentum. This leads to transition 
from WL at low density to WAL at higher densities in 2D 
hole systems caused by increase of spin-orbit strength.^! 

Experimental studies showed the WL-WAL transition 
with increase of 2D hole density in both temperature and 
magnetic-held dependences of the conductivity.®^ The 
WL theory for hole systems^ffi^ was developed for the 
so-called "diffusion" regime where the backscattering oc- 
curs after many collisions with impurities. The diffusion 
theory describes well the magneto conductivity in low- 
mobility hole structures,^ but it is inapplicable to high- 
mobility systems studied nowadays because it describes 
magnetoconductivity in very low helds. 

The temperature dependence of the conductivity is 
also different in high-mobility systems. WL is present due 
to these closed paths where the interference is not bro- 
ken by dephasing processes. This means that the phase 



breaking time should be much longer than the momen- 
tum scattering time r. However, the diffusion theory of 
WL assumes a stronger condition In (t^/t) 3> 1 which is 
not realized in high-mobility structures with long scat- 
tering times. The backscattering in these systems occurs 
after carrier propagation along "ballistic" paths with few 
impurities which are ignored by the diffusion theory. The 
aim of this work is to develop the WL theory of two- 
dimensional holes valid in both the diffusion and ballistic 
regimes and calculate the temperature dependence of the 
conductivity correction. 

The hole states in quantum-well structures grown from 
A3B5 materials are described by the Luttinger Hamilto- 
nian where the wavevector component along the growth 
direction is an operator. We apply the isotropic approx- 
imation and assume the quantum well to be rectangular 
and infinitely deep. We have found the hole energy dis- 
persion in the subbands of size quantization in this model 
by a standard wayP^ The dependence of the hole energy 
in the ground subband on the two-dimensional wavevec- 
tor, E(k), is non-parabolic at ka > 2, see inset to Fig.JT] 

The size-quantized hole levels in a symmetrical quan- 
tum well are double degenerate at each value of k. We 
choose these two states as symmetric (s) and asym- 
metric (a) relative to a mirror reflection in the sym- 
metry plane of the quantum wellP^ The corresponding 
wave functions are ^! a k = F ak {z) exp (ife • p), where z 
is the direction of quantization, p is the in-plane co- 
ordinate, a — s, a, and F ak {z) are linear superposi- 
tions of the four Bloch functions u m of the total angu- 
lar momentum 3/2 and projections on the growth axis 
m = ±3/2, ±1/2: u ±3/2 = =f(A ± iY)s±/\/2, u ±l/2 = 

[T(X ± iY)s T /V2 + V2Zs ± ]/V3. Here X, Y, Z are the 
Bloch functions of the top of the valence band, and s± 
are spin functions with projections ±1/2 on the z axis. 
In the basi s \u3j 2, U1/2, U-1/2, U-3/2], F ak have the fol- 
lowing for m 1 11 ! 16 -! 

F sk (z) = (1) 
[~voC(z) 1 w 1 S{z)e i ^,-v 2 C(z)e 2i ^,w 3 S(z)e 3i ^] , 
F ak {z) = (2) 
[iv 3 S(z)e' 3i ^ ,v 2 C(z)e- 2i ^ ,i Vl S{z)e^ ,v Q C(z)] , 
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FIG. 1: Elastic relaxation times of holes. Inset: energy dis- 
persion of the ground hole subband in a 100 A wide quantum 
well. 



where ipj, is the polar angle of the wavevector fc. The 
dependence on fc of the real coefficients K0...3 as well as 
symmetric (C) and antisymmetric (S) functions of the 
coordinate z is determined by the energy dispersion E{k). 

We consider elastic scattering of two-dimensional holes 
by the short-range potential V{r) — Vo^ Ji 8{r — ri), 
where r$ = (p^, z{) are coordinates of the impurities. The 
spin-orbit coupling leads to anisotropy of scattering even 
by such isotropic impurities. Indeed, the hole wavefunc- 
tions Eqs. ([l}, ^ depend on a direction of the wavevector 
fc. As a result, the matrix elements of scattering from the 
state \ak) to the state \j3k') by the potential V(r) 



(h-h')-pi 



depend on the angular coordinates of the initial and final 
wavevectors ifk and ifk 1 ■ Therefore the elastic scattering 
probability contains the 1st, 2nd and 3rd Fourier har- 
monics of the scattering angle 9 = ifk> — y>k in addition 
to the isotropic term. Hence the scattering times de- 
scribing decay of the nth Fourier harmonics of the hole 
distribution function 

T n h J Ztt 


are different for n = 1,2,3. All higher harmonics decay 
with the time tq given by 
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2tt r de 



9 / 7j-(\v..(0)\ 2 + \Va.(e)\ 2 ) 



h J 2tt 




Here angular brackets denote averaging over the impu- 
rity positions, and g = k/(2nhvk) is the density of states 
where Vk = h~ 1 dE/dk is the velocity of holes with the 



wavevector fc. The dependences r„(fc) for n = 0,1,2,3 
are plotted in Fig. [I] for homogeneous impurity distri- 
bution over the z coordinate (averaging is restricted to 
integration over Zi). Figure [T] demonstrates that the 
p-type heterostructures differ qualitatively from topo- 
logical insulators and graphene where scattering is also 
anisotropic, but the ratios T\/tq — 2, t„/t = 1 (for 
n > 1) are fixed, and all relaxation times are indepen- 
dent of fc. In contrast, the spin-orbit coupling is ab- 
sent for holes at fc — > 0, where the spectrum is parabolic 
and scattering is isotropic, while at ka > 2 both condi- 
tions are violated. This implies the spin-orbit coupling 
strength in the hole systems is controlled by the density 
p which sets the Fermi wavevector fcp = y/2wp. 

This anisotropy of scattering complicates the calcula- 
tion of the WL conductivity correction. It is determined 
by the Cooperon C"g (<Pk><Ph' ,q) which depends on the 
four indexes enumerating the hole states (a, /3, 7, 8 = 
s, a), on the directions of the momenta at the Fermi circle 
VfciVfc'i an d on the two-dimensional vec tor q, q <C fcp. 
The Cooperon satisfies the equatio n 1 11 ! 17 ! 
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where 



P(<Pi,q) = 



1 



2TTgpTo I H 
iqv F T cos (ipi - tp q ) '' 



gp and vp are the density of states and the velocity at 
fc = fcp, and the correlators are defined as 

For scattering in a symmetrical quantum well the corre- 
lators have the properties following from Eqs. ([I]), ([2]): 
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and all other correlators are equal to zero. Note that the 
relaxation times are determined by (\V SS \ 2 + |V as | 2 ) — 
~ ^as- The correlators contain Fourier har- 
monics exp (±imy>k ± inipk' ) with m,n = 0...6, and 
Wgs > W B a ) depend on the difference cpk ~ <Pk' while 
Wsaifk, fk') does not. 

The correlator properties Eqs. (|5) allowed us to de- 
couple partly Eqs. |3| and derive independent equations 
for the Cooperons C + and C_ as well as two coupled 
equations for the Cooperons C\ and C2 defined as 
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The equations are as follows: 



(7) 



271 



2ir 



P(ip u q) TO((9i) T WTO)] C ± (tp u ipv, q), 



C 1 (<p h ,<p v ,q) = W^{0) 



+W s s Z(vk,Vi)C 2 (<p 1 ,<p k >,q) 



(8) 




FIG. 2: Conductivity correction as a function of hole density 
at t(O)/t-0 = 0.01,0.02,0.04 and 0.09. Inset: the backscatter- 
ing (a a ) and non-backscattering (ctj) contributions. 



(9) 



Here we introduced 9 = ip^ — tp^ and 9\ = p^ — Lp\. 

We have found the Cooperons C+, C_, Ci and C*2 
expanding them in Fourier series which transfers the in- 
tegral Eqs. ([7|-([9]) to systems of linear equations. Since 
the correlators have a finite number of Fourier harmon- 
ics, the linear equation systems are finite: 7 independent 
equations for both C+ and C_, and 14 coupled equa- 
tions for C\ and C%. The Cooperon Fourier harmonics 
have been found by a numerical solution of these linear 
equation systems. 

The found Cooperons allowed us to calculate the WL 
conductivity correction a. It is well known that this cor- 
rection is given by a sum of two maximally crossed di- 
agrams. Accordingly, it equals to a sum of two terms, 
(7 = a a + er&, where a a arises from the backscattering 
processes and tr& contribution is due to coherent scatter- 
ing by arbitrary angles. These corrections are expressed 
via the Cooperons as follows 
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Here 9 = ip — ip', the upper index (2) and (3) indicates 
that the diagrams are started from 2 and 3 impurity lines, 
and the factor r 2 comes from the renormalization of two 
vertexes due to the scattering anisotropy. 

The results of calculation of the WL conductivity cor- 
rection is shown in Fig. [2] The WL to WAL transi- 
tion with increase of the hole density is clearly seen for 
t(O)/t0 < 0.04. It occurs at kpa « 2 which corresponds 
to the density p s» 7 x 10 11 cm -2 for the quantum well 
width a = 100 A. However, Fig. [| demonstrates that the 
conductivity is a non-monotonous function of the hole 
density which changes sign and become again negative at 
high k^a for shorter dephasing times r^. Inset in Fig. [2] 
shows the backscattering and nonbackscattering contri- 
butions, a a and 

The backscattering contribution dominates for all val- 
ues of the parameter kpa. This can be explained by the 
approximate relation 



T\ - r 

( 

n 



which is derived in the diffusion approximation.^^ The 
difference between t\ and To is nonzero, but it does not 
exceed 30%, see Fig. [I] Therefore the nonbackscat- 
tering contribution crj does not reach the large value 
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FIG. 3: The backscattering correction at t(0)/t^ = 0.07 cal- 
culated exactly (solid line) and in the diffusion approximation 
(dashed line). Inset: spin relaxation rates. 
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FIG. 4: The temperature dependence of the conductivity cor- 
rection at different hole densities. Inset highlights the region 
where the conductivity has a minimum. 



(e 2 /2Tr 2 h) In (t^/ti). In contrast, it does not exceed the 
value at k = where o\, = (e 2 /2ir 2 h) In 2, see Fig. [2] 

In order to explain qualitatively the non-monotonous 
behaviour of the conductivity correction with the hole 
density we also compare the results of our exact calcula- 
tion with ones obtained in the diffusion approximation: 
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Here the longitudinal and the transverse spin relaxation 
rates are given by^ 
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The dependence Eq. ( 12 ) is shown in Fig. [3] together with 
the result of the exact calculation. The increase of the 
WL correction at < kpa < 3 is described by fast in- 
crease of both spin relaxation rates which makes smaller 
two first logarithmic terms in Eq. (12). This demon- 



strates that WL to WAL transition in hole systems at 
long dephasing times is driven by increase of spin-orbit 
coupling with increase of the hole density. However, 
at k^a > 3 the WL correction calculated in the diffu- 
sion approximation saturates at the value smaller than 
(e 2 /Ait 2 K) In (t^/ti) obtained in the limit of very fast spin 
relaxation. This is caused by decrease of the transverse 
spin relaxation rate 1/t s j_ and by saturation of the lon- 
gitudinal rate l/r s ||, see inset to Fig.[3) 

This behaviour of the spin relaxation rates at ka > 
1 can be explained by considering limit of a very wide 



quantum well. The coefficients Wo.. .3 m Eqs. Q, Q in 
this limit are i>o = 1/2, V\ = v% — 0, i>2 = v3/2p^ so the 
symmetrical and asymmetrical wave functions have the 
following forms: 



nz\ exp (±iipi~) 



[(Fcos (fk - Xsmipk)s± + iZe ±iipk s T ] 



These wave functions correspond to the angular mo- 
mentum projections ±3/2 on the direction of k in the 
structure plane as it should be at large ka. Calculating 
the scattering matrix elements V ss — Vo(F s k'\F s k) and 
V sa = Vo(F s k'\F a k) we derive that at ka — > 00 

1/t s± = 0, 0.9r s || =t = t 1 . 

Infinite t s j_ shows that the hole spin relaxation is sup- 
pressed not only at small but also at large ka. According 
to the calculation results shown in inset to Fig. [3] this 
limiting case is almost realized at k-pa < 4 when only the 
ground subband is occupied. Therefore the first logarith- 
mic term in Eq. ( 12 1 is small at kpa ~ 4 while the second 



one increases partially compensating the last term, and 
the total correction saturates. Moreover, in contrast to 
the diffusion-approximation result, the exactly calculated 
WL correction totally repeats the behaviour of ti/t s ±: 
the above derived suppression of hole spin relaxation in 
wide quantum wells is reproduced in the WL conductiv- 
ity correction as it is seen from exact calculation result 
shown Fig. [3] which has a maximum and decreases at 
larger ka. 

We have also calculated the WL conductivity correc- 
tion as a function of temperature T ~ tT 1 at different 
hole densities. The result shown in Fig. [4] demonstrates 
the transition from metal to insulator temperature be- 
haviour. This transition is due to spin-orbit interaction 
which is present even in symmetrical quantum wells with 
the degenerate hole energy spectrum. Inset to Fig. [4] 
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shows non-monotonous temperature dependence at mod- 
erate hole densities. 

To summarize, we have developed WL theory for holes 
in quantum wells taking into account both ballistic and 
diffusion processes. The conductivity correction is cal- 
culated in a wide range of hole densities and tempera- 
tures taking into account real nonparabolicity of the en- 
ergy spectrum. WL instead of WAL can occur at high 
densities due to suppression of hole spin relaxation at 
large in-plane wavevectors. The transition from metal 
to insulating temperature behaviour is demonstrated in 
the symmetrical quantum well with spin-degenerate hole 



subbands. Such transition can be detected by magne- 
toconductivity measurements. This mechanism can con- 
cur with the mechanism due to Rashba splitting of hole 
energy spectrum which has been used in Ref. 121 for ex- 
planation of WL experiments on two-dimensional hole 
systems. 
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